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Abstract 

We investigate superfluid properties of a one-component Fermi gas with a uniaxially anisotropic 
p-wave pairing interaction, Ux > Uy = Uz (where Ui {i = x, y, z) is a pj-wave pairing interac¬ 
tion). This type of interaction is considered to be realized in a Fermi gas. Including pairing 
fluctuations within a strong-coupling T-matrix theory, we determine the px-wme superfluid phase 
transition temperature TA, as well as the other phase transition temperature (< TA), 

below which the superfluid order parameter has the px + ipy-wave symmetry. In the normal state 
near Tf "", px-wave pairing fluctuations are shown to induce an anisotropic pseudogap phenomenon, 
where a dip structure in the angle-resolved density of states around a; = 0 is the most remarkable 
in the Px direction. In the Px-wave superfluid phase < T < TA), while the pseudogap in 

the Px direction continuously changes to the superfluid gap, the pseudogap in the perpendicular 
direction to the px axis is found to continue developing, because of enhanced py-wave and pz-wave 
pairing fluctuations around the node of the px-wave superfluid order parameter. Since pairing 
fluctuations are always suppressed in the isotropic s-wave superfluid state, this phenomenon is 
peculiar to an unconventional Fermi superfluid with a nodal superfluid order parameter. Since the 
p-wave Fermi superfluid is the most promising non s-wave pairing state in an ultracold Fermi gas, 
our results would contribute to understanding how the anisotropic pairing fluctuations, as well as 
the existence of plural superfluid phases, affect many-body properties of this unconventional Fermi 
superfluid. 

PACS numbers: 03.75.Ss,05.30.Fk,67.85.-d 
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I. INTRODUCTION 


An interesting feature of a non s-wave Fermi condensate is that it may have plural 
superfluid phases, originating from active orbital and/or spin degrees of freedom. Indeed, 
this possibility has experimentally been conhrmed in various Fermi superfluid systems, such 

[l|, as well as superfluid liquid ^He j^, y] . However, 
lat has already been realized is still only the simplest s-wave one in 
. Thus, going beyond this situation is a crucial challenge in this 


as heavy fermion superconductor UPts 
the pairing symmetry t 
cold Fermi gas physics 
research held. 


In this regard, a p-wave superhuid Fermi gas is a strong candidate, and the possibi 


of this spin-triplet 
and theoretically 


miring state has extensively been discussed both experimentally js 


ity 


resonance in 


m 


19M3 

lIllQ 


33| . Several experimental groups have discovered a p-wave Feshbach 
10] and ®Li [ll, 12] Fermi gases, so that we can now tune the strength 


of a p-wave interaction from the weak-coupling regime to the strong-coupling regime, by 
adjusting an external mag netic held. This experimental development has realized p-wave 


Feshbach molecules 


11 


14l-ll7l|. Thus, although one still needs to overcome some diffi¬ 
culties, such as the three-body loss j^, 34, 3) dipolar relaxation InlQ 

(that destroy p-wave molecules [3^), the p-wave superhuid state seems a very promising non 
s-wave pairing state in an ultracold Fermi gas. 

It has been predicted 2, 27] that, when a p-wave interaction has a uniaxial anisotropy, 
a p-wave superhuid Fermi gas may have two superhuid phases with diherent p-wave pairing 
symmetries. Such an anisotropic p-wave pairing interaction is considered to be realized in a 


Fermi gas , because the split of a p-wave Feshbach resonance into a p^-wave channel 
and degenerate Py-wave and p^-wave channels by a magnetic dipole-dipole interaction has 
been observed, when an external magnetic held is applied in the x direction. Since the 
observed resonance held of a p^-wave Feshbach resonance is higher than that of the other 
degenerate channels 10| , a p-wave pairing interaction associated with this p-wave Feshbach 
resonance has the uniaxial anisotropy, that is, a pa;-wave pairing interaction Ux is stronger 
than py-wave {Uy) and p^-wave (t4) interactions. As a result, the Pa,-wave superhuid state 
has the highest superhuid phase transition temperature In this case. Refs. Q, 


27 


have 


pointed out that the system experiences the other phase transition from the Px-wave state 


to the Pa; -t-ipy-wave state at Tc 


,Px+iPy 


< when the uniaxial anisotropy satishes a certain 
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condition. Thus, the realization of a p-wave superfluid Fermi gas would enable us to study 
physics of plural superfluid phases, from the weak-coupling regime to the strong-coupling 
limit in a systematic manner. 

When the above-mentioned p-wave superfluid Fermi gas is realized, the existence of strong 
pairing fluctuations near jg an interesting research topic. In the Px-wave superfluid 

phase, since single-particle excitations are gapless in the nodal direction (T px) of the px- 
wave superfluid order parameter Ap^{p) oc Px, Py-wave and p^-wave pairing fluctuations can 
continue developing, to be the strongest at _ Thus, even far below TJ’"', an anisotropic 

pseudogap phenomenon is expected in the nodal direction of the Px-wave superfluid order 
parameter near pP^'^^Py _ 

In an isotropic s-wave superfluid Fermi gas, since the BCS gap opens in all the momentum 
direction of single-particle excitations, pairing fluctuations are soon suppressed below the 
superfluid phase transition ternperature T^.. Indeed, it has been shown that a pseudogap in 
the density of states above Tc |37 h 44| soon changes to the s-wave BCS superfluid gap below 
Tc 451]. Even in the px-wave case, the enhancement of pairing fluctuations around the node 
would not occur in the absence of a Py-wave and a px-wave interactions, because px-wave 
pairing fluctuations are soon suppressed by the Px-wave superfluid order below Thus, 
the above-mentioned pseudogap phenomenon in the px-'wave state near j'P^^'^Py jg peculiar 
to an unconventional Fermi superfluid with a nodal superfluid order parameter, as well as 
with plural superfluid phases. 

In this paper, we theoretically investigate strong-coupling properties of a one compo¬ 
nent p-wave superfluid Fermi gas with a uniaxially anisotropic p-wave pairing interaction 
[Ux > Uy = Uz). Including p-wave pairing fluctuations within the framework of a strong 
coupling T-matrix approximation, we determine Tp^ and pP^'^^Py _ in the p^-wave superfluid 
phase, we calculate the angle-resolved single-particle density of states, to clarify that pairing 
fluctuations in the nodal direction (T px) of the px-'wave superfluid order parameter con¬ 
tinue to develop, leading to a pseudogapped single-particle excitation spectrum in the nodal 
direction. 

The paper is organized as follows. In Sec.II, we explain our strong-coupling formalism 
for a p-wave superfluid Fermi gas. In Sec.Ill, we show our numerical results on Tp^ and 
pP^'^^Py _ Here, we clarify the condition for the appearance of Px-wave and Px + ipy-wave 
superfluid phases. In Sec.IV, we examine the angle-resolved single-particle density of states 
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(ARDOS). We show that a pseudogap anisotropically appears in ARDOS, not only in the 
normal state near but also in the px-weive superfluid phase near Xo characterize 

this anisotropic many-body phenomenon, we introduce the characteristic temperature T* as 
the temperature below which a dip structure appears in ARDOS. Throughout this paper, 
we take h = ks = 1, and the volume of the system V is taken to be unity, for simplicity. 


II. FORMULATION 


We consider a one-component Fermi gas with a uniaxially anisotropic p-wave interaction, 
described by the Hamiltonian, 


H 'CpCpCp 2 ^n{p)PiUiP^Fn{p)c^p_^_^^2^^-p+q/2^-p'+Q/‘^^P'+Q/'^- (^) 

p p,p',q i=x,y,z 

Here, is the creation operator of a Fermi atom with the kinetic energy = Sp — p = 
p^/2m — /i, measured from the Fermi chemical potential p (where m is the atomic mass). 
—Ui (< 0) is a pairing interaction in the pj-wave Cooper channel {i = x,y,z), having the 
uniaxial anisotropy Ux > Uy = [l0|. In this paper, we do not deal with details of a 

p-wave Feshbach resonance, but simply treat {Ux,Uy,Uz) as a tunable parameter set. To 


eliminate t 


2G-22 


24 


le ultraviolet divergence, the last term in Eq. ([T]) involves the cutoff function 

m, 


Fn{p) = 


( 2 ) 


1 + (P/Pc)^"’ 

where Pc is a cutoff momentum. For simplicity, we use the same cutoff function Fn( p) f or all 
the p-wave Cooper channels. Equation (12|) gives a Lorentzian cutoff when n = 1 241) 


gives a sharp cutoff when n = oo 
in Sec. HI. 

As usual 


33| . We will discuss the cutoff dependence of our results 


20 


22 


24 


33| , we conveniently measure the strength of a p-wave interaction 


m 


terms of the p^-wave scattering volume Vi {i = x, p, z), as well as the effective range ko, that 
are given by, respectively, 

Anvi Ui 1 


m 


U, 


E 


„2 : 


p t' 


( 3 ) 


( 4 ) 
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Since we take the same cutoff function Fn{p) for all the pj-wave interaction channels, the 
effective range fco are channel-independent. We also introduce the anisotropy parameter, 




W (> 0)- 


(5) 


Then, the p-wave interaction can be specihed by the parameter set {v~^,6v~^,ko,n). The 
weak-coupling side and the strong-coupling side are characterized as ^ 0 and 

(PfW)~^ ^ 0, respectively, where pp is the Fermi momentum. 

To deal with strong-coupling phenomena in the superfluid phase, it is convenient to 
rewrite the model Hamiltonian in Eq. (Ip) i nto the sum of the mean-held BCS part and the 


term describing huctuation corrections 
have. 


45 


46[ |. Under the Nambu representation we 


where 




( 6 ) 


q,t=x,y,z 




(7) 




is the two-component Nambu held, and r* {i = x, y, z) are Pauli matrices acting on particle- 


hole space 


47l |. The hrst term in Eq. ([6]) is just the mean-held BCS Hamiltonian, where 


A(p) = 


0 A(p) 
A*(p) 0 

is a 2 X 2 matrix p-wave superhuid order parameter. Here, 


( 8 ) 


A(p) = b-pFn{p), 

where b = (bx^^by^bz) has the form, 

bi Uj ^ (p) (C—pCp) . 


(9) 


( 10 ) 


The second term in Eq. ([6]) gives huctuation corrections to the mean-held Hamiltonian, 
where the so-called generalized density operator [^, [4^ , 


pfiQ) = , 


( 11 ) 


5 
















(a) 2 ip, /co„) 



FIG. 1: (a) Self-energy correction 'S{p,iujn) in the T-matrix approximation (TMA). The wavy line 
is the particle-particle scattering matrix [q^ivn) given in (b). The solid line and the dashed 
line describe the mean field single-particle Green’s function Go{p, ioJn), and the p-wave interaction, 
respectively. In (a), the factor {pi — qi/2)Fn{p — is assigned to each vertex (solid circle). In 

(b), we assign piFn{p)Ts to each vertex. In this figure, —s means the opposite sign to s = ±. 


physically describes pi-wave superfluid fluctuations (where r± = (ri ± iT 2 ) /2). 

Strong-coupling corrections to single-particle excitations can be conveniently described 
by the self-energy S(p, in the single-particle thermal Green’s function, 

1 


G{pi iuiYn) 


^0 ^ iPi 


( 12 ) 


where is the fermion Matsubara frequency. In Eq. flT^ . 

1 


Go{p,i(^m) = - 


il^m - ^pTs + A(p) 


(13) 


level. Treating the last term in 


is the single-particle Green’s function in the mean-£e 
Eq. ([6]) within the T-matrix approximation (TMA) 40 
expression for the self-energy S(p, ia;„) shown in Fig. [T](a) |48|, which gives 


41 


4^, we have the diagrammatic 


i,j=x,y,z q,iUn 



r Ti 

r m 

2/ 

r “ 

b - 2 J 


X 


GoiP - Q, - Wn)ri/(q, iVn) Gf{p - q, - w„)r++(qr, iUn) 

Gl'^ip - q, iUJn - (g, '^^n) 


( 14 ) 
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Here, Vn is the boson Matsubara frequency. The TMA particle-particle scattering matrix 
(qr, iun) in Eq. ffTTl) obeys the equation, 


Ti’j{q,iUn) = -Ui6ij6s,-s'-U^Y^ n";f (15) 

s"=± k=x,y,z 

where —s means the opposite sign to s = ±, and 


n-;/(g, ii^n) = ^ ^PiPjF,?(p)Tr TsGq (p + iujn) vGo (p - 


Q 


q . 


I Ur 


(16) 


is the pair correlation function. In particular, H^ {q, iun) describes fluctuations in the 
Pi-wave Cooper channel, and {q,Wn) (i 7^ j) describes coupling between pj-wave and 
Pj-wave pairing fluctuations. 


In the present anisotropic case (t4 > Uy = Uz), the highest superfluid phase transi 


temperature is obtained in the Pa;-wave Cooper channel. As in the s-wave case 


45 


ion 


m 


the TMA gap equation for the p 3 ;-wave superfluid order parameter Ap^(p) = h^PxFnip) is 
obtained from the Thouless criterion [rP*),,(qf = 0,zz/„ = 0)] ^ = 0 in the Pa;-wave Cooper 
channel, where 


= T + r.+ (q, iv„) - p, {q,iv„) - r++(ij, !!/„)] , 


(17) 


describes the phase fluctuations of the Pa;-wave superfluid order parameter Ap^(p). Physi¬ 
cally, this guarantees the existence of a gapless Goldstone mode associated with the broken 
U{1) gauge symmetry. Noting that 6^ = 6^ = 0 in the p^-wave superfluid phase, we obtain 
the gap equation in the Pa;-wave superfluid state as 


127 rna, 2 172/ \ 

1 = - 

p 


m 


: tanh 


^p+ W.ipW 

2T 


2^ei+wjp)\- 

The equation for is obtained from Eq. (1T8|1 by setting Ap^ (p) = 0 as 

^ ■ tanh ^ 


2er 


(18) 


127rr’3, 2 172/ \ 

1 = — Yp^^nip) 

p 


m 




2Tt 


2er 


(19) 


In the case of uniaxially anisotropic p-wave interaction. Ref. 27| pointed out that, without 
loss of generality, one may restrict the structure of the p-wave superfluid order parameter 
to the form. 


( 20 ) 



b. 


Bx 

b = 

by 

= 

iBy 


1 i 


0 , 
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where and By are real quantities. Thus, in the px-wecve superfluid phase below (where 
Bx = bx 0 and By = 0), the other possible superfluid instability is only associated with 
the Px + ipy-weive one, having the superfluid order parameter, 

^Px+ipy{p) ~ [BxPx 'i'ByPy]Fn{p). (21) 


Since Bx is already present below the superfluid phase transition temperature 

rjnPx+tpy jg determined from the Thouless criterion 491 in the Py-wscve Cooper channel 

[rph(q = o,*z/„ = o)]"' = o. 


1 = 


m 


: tanh 


Q + WxipW 1 

2^rjnPx+ipy 


( 22 ) 


where the pa;-wave superfluid order parameter Ap^(p) = bxPxFnip) obeys the gap equation 

ca. 

We numerically solve Eqs. (ITSll . ([T9]), and (12^ . to self-consistently determine 
and Ap^(p). In this procedure, we also solve the equation for the total number Ap of Fermi 
atoms. 


Nv = 


T 


F, Tr|TsG(p,!w„)], 


(23) 


p,tU}n 


to include strong-coupling corrections to the Fermi chemical potential p. 

We examine the anisotropic pseudogap phenomenon by calculating the angle-resolved 
single-particle density of states (ARDOS), 


^OO 9 


P{.^,P) = — 


TT 


p^dp 

(27r)3 


Im [G 11 (p, icUn —)■ CJ -h fh)] , 


(24) 


where p = p/|p|, and Gii{p,i(^n oJ + iS) is the (1,1) component of the analytic continued 
TMA Green’s function in Fq. flT^ . ARDOS in Eq. fl2T|) is related to the ordinary density 
of states p{u) as 

p{u) = j sin 9pd9pd(l)pp{u,p). (25) 

Here, we choose the px axis as the polar axis {px = p cos 9p). 

We note that, since the anisotropy of the px-wave superfluid order parameter Ap^(p) cx Px 
lowers the symmetry of the system, we need much time to compute the number equation 
fl2^ below compared to the symmetric s-wave case. To avoid this difficulty, in this 
paper, we approximate the TMA Green’s function G{p, iujn) in the number equation (|23|) to 


G(p, iUn) ^ Go(p, iuJn) + Gq{p, iUn)t{p, ia;„)Go(p, iu:n)- 


( 26 ) 











Equation fl26|) is just the same form as the Green’s function in the strong-coupling theory 
developed by Nozieres and Schmitt-Rink (NSR) The NSR theory has extensively been 
used in the s-wave case, to successfully explain the BCS-BEC crossover behavior of the 


superfluid phase transition temperature 
in the crossover region 


50 


with Ux = Uy = Uz 


52 


51| . as well as the superfluid order parameter 


53| . The NSR theory has also been extended to the p-wave case 
20|]. Thus, we expect that the NSR Green’s function in Eq. (|26|) also 


works in determining , /i and Ap^(p). 

On the other hand, it is also known that the NSR theory unphysically gives negative 


density of states in the BGS-BEG crossover region 


41 


54| . Since this serious problem is 


absent in TMA, we use the TMA Green’s function in Eq. fll2p . in considering single-particle 
properties of a Px-waye Fermi superfluid. 

Here, we summarize our detailed parameter settings. For the effective range, we take 
ko = —SOpF, following the experimental result on a Fermi gas |^. For the cutoff 
function Fn{p) in Eq. ([2]), we set n = 3. The cutoff momentum pc in Fn{p) is determined 
so as to reproduce ko = —SOpp, which gives Pc = 27pF. Since we only deal with the normal 
state, as well as the px-wave superfluid state, ARDOS in Eq. (j2^ is actually independent 
of the angle (j)p around the Px axis. Thus, the anisotropy can be simply specihed by the 
polar angle 6p measured from the px axis. Noting this, we write Eq. fl2T|) as p{u, Op) in what 
follows. 


III. PHASE DIAGRAM OF AN ULTRACOLD FERMI GAS WITH p-WAVE IN¬ 
TERACTION 


Figure [2] shows the phase diagram of a one component ultracold Fermi gas in terms of the 
p-wave interaction strength, {ppVx)~^, and the temperature. When 6v~^ = 0 (f4 = Uy = 
Uz), Fig. [2](a) shows that the superfluid phase is dominated by the px + ipy-wave pairing 
state. This superfluid region gradually shrinks with increasing the uniaxial anisotropy 
as shown in Figs. [2](b) and (c). Since the present anisotropy {Ux > Uy = Uz) favors the 
Px-wave symmetry, the region of the px + ipy-wave state eventually vanishes, as shown in 
Fig. El^d). We briefly note that the overall structure of this phase diagram is consistent with 
the previous work based on mean-held analyses 0,0- 

In addition to the coupled equations (IT^ and fl23|) also give the Fermi chemical 


9 
























FIG. 2; (Color online) Phase diagram of a one-component Fermi gas with a uniaxially anisotropic 
p-wave pairing interaction. The solid line and the dashed line are and ^ respectively. 

T-p is the Fermi temperature. 


potential ix{T = Tf) shown in Fi^ |3l which exhibits the typical BCS-BEC crossover 


behavior 


20 


21 


33 


3 , 0 , 


45 


50 


52]. That is, with increasing the interaction strength, 


p(r = gradually deviates from the Fermi energy Sp, to be negative in the strong¬ 
coupling regime, when {ppVx)~^ ^ 0 . 

At .j-pg coupled equations flT 8 |) . fl22|) and (l2^ also give /i(T = as 

well as the pa;-wave superfluid order parameter Ap^{p,T = 7 "^^+*^!'). Figure [3] shows that 
fi{T = ~ ^(T = in the whole interaction regime, indicating that the chem¬ 

ical potential is almost T-independent in the px-'weive superfluid phase. For the pa;-wave 
superfluid order parameter Ap^{p,T = 7 '^“'+*^*')^ of course, this quantity has already existed 
above as shown in Fig. 01 Although the pairing symmetry changes from the px-wave 

one to the Px + ipy-wave one at 7 ^^+®^’*'^ it is known 27 1 that this symmetry change occurs 
smoothly, in the sense that By in Eq. fl^Ill continuously grows from zero below 7 ?^+*^’*'. 
Thus, the second order phase transition is expected at 7 ^^+®^*' (unless the superfluid order 
parameter exhibits an unexpected discontinuity at 7 ^“'+®^*'). 

In Fig. m we see that Ap^{p) has a discontinuity at which is, however, an artifact 
of TMA we are using in this paper. The same problem has already been known in the 


s-wave case 


45 


53 


55| . In the latter case, it has been pointed out 


55 II that one needs to 


correctly include an effective repulsive interaction between Cooper pairs beyond TMA, in 
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FIG. 3: (Color online) Calculated Fermi chemical potential /x at the two superfluid phase transition 
temperatures, Ti"" and . We set (p|(5u)“^ = 0.3. 



0.02 0.04 0.06 0.08 0.1 


r/7p 


FIG. 4: (Color online) Calculated factor bx in the px-wave superfluid order parameter Ap^ = 
bxPxF'n= 3 {p), when {pp6v)~^ = 0.3. Each result ends at _ 

order to recover the expected second order phase transition. Although this improvement 
is also crucial in the p-wave case, we leave this problem as a future problem, and examine 
strong-coupling effects in the p 3 ;-wave superfluid phase within TMA. 

In the weak-coupling BCS limit, the number equation fl25]) simply gives p = ep, so that 
the superfluid phase transition temperature is determined from Eq. flT^ with p = e-p- 
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,(7Vf/2) 


,(AAf/6) 


(pJSu)-! 


FIG. 5: (Color online) Calculated two superfluid phase transition temperatures, and Tp^^ip^, 
as functions of the anisotropy parameter {pp6v)~^. In panel (c), the solid circles shows the BEC 
phase transition temperature Tbec obtained from Eqs. ([2711 - (l^^ . The upper dashed line in this 
panel shows Tbec(-^f/ 2). The lower dashed line shows Tbec(-^f/6). 


The resulting does not depend on the anisotropy parameter for a hxed pa;-wave 

interaction strength Vx- On the other hand, gradually comes to depend on as one 
approaches the strong-coupling regime, as shown in Fig. |5l 

To explain the anisotropy dependence of in the strong-coupling regime shown in Fig. 
|5](c), we first note that the system in the strong-coupling limit ^ may be viewed as an ideal 
gas mixture of three kinds of tightly-bound molecules (with the molecular mass M = 2m) 
that are formed by the three pi-wave interactions {i = x, p, z). Indeed, in the strong-coupling 
limit, the number equation fl2^ at is reduced to 


Nf 




(27) 


where 

^b = E’‘b(^-Fb) (28) 

<7 ^ 2 

is the number of molecules in the Pi-weive Cooper channel, with nB((n) being the Bose 
distribution function. The TP'^-equation flT^ gives the Bose chemical potential pg in Eq. 


(p% = 0, 

I Yb = Yb 


26v ^ 




(< 0). 


(29) 
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In the absence of uniaxial anisotropy ((5n“^=0), all the three components simultaneously 
satisfy the BEC condition, /^b = hn = /^b = 0- Thus, the phase transition temperature 
(= Tbec(^f/6 )) is determined from the equation, iVp/b = = iVg = iVg, which gives 

‘2,71 f N \ 

rBEc(iVF/6) = {~f) “ 


where C(3/2) = 2.612 is the zeta function. 

In contrast, when > 0, Eq. fl2^ shows that the Py- and pz-wave components no 
longer satisfy the BEC condition. In the extreme case when dv~^ 3> 1, the Bose chemical 
potentials /ig and /ig in Eq. fl29l) are much lower than zero, so that one can ignore the 
contributions of these components to the number equation flTTI) . as Np/2 = N§. This 
means that most atoms form bound molecules in the Pa,-wave Cooper channel, which is 
quite different from the case of = 0, where only one third of Fermi atoms contribute to 
Pa;-wave molecules. Because of this, the BEC phase transition temperature (= Tbec(-^f/ 2)) 
in this extreme case is higher than of = 0 in Eq. fl30|) as, 

Tbeo(A'f/ 2) = =0.137rF. (31) 


Figure |5]^c) shows that the BEC phase transition temperature Tbec(-^b) calculated from 
Eqs. fl?n) - fl2^ monotonically increases from Tbec(-^f/6 ) to Tbec(-^f/ 2), with increasing 
the anisotropy parameter The well agreement of with Tbec shown in this hgure 

indicates that the anisotropy dependence of in this regime comes from the increase of 
the molecular bosons in the Px-weive Cooper channel, with increasing the uniaxial anisotropy 
of the p-wave interaction. 

In contrast to we see in Fig. [5] that the px + ipy-wave superfluid phase transition 
temperature decreases with increasing 6v~^, to eventually vanish at a critical value 

(Note that this vanishing has already been expected in Fig. [2j) Evaluating 

this critical value in the whole interaction strength, we obtain Fig. |6l This hgure shows 
that is not so sensitive to the interaction strength, to always lie in the narrow range, 
0.4 < (p|5n)“^ < 0.6. 

To understand this behavior of since thermal huctuations are absent at T = 0, it is 


convenient to employ the BCS-Leggett theory 


57l |. which consists of the coupled Eq. 
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FIG. 6: (Color online) Critical value Sv~^ of the anisotropy parameter at which vanishes 

(solid line). The dashed line shows the result within the BCS-Leggett theory. In obtaining the solid 
line, we have taken a small but finite value of T ( < O.OITp) because of computational problems. 


at = 0 with the mean-field number equation at T = 0, 


Y Tr[r3Go(PC^^n)] 

P,i^n 



1 - 




(32) 


As shown in Fig. El the BCS-Leggett theory semi-quantitatively reproduces the TMA result 
for weak-coupling BCS regime (|Ap^(p)| <C ep), the number equation fl5^ 

simply gives fi = £p. Substituting this into Eq. fl2^ with = 0, one obtains the upper 

bound of the BCS-Leggett theory as 


= - = 0.64. (33) 

TT 

In the strong coupling regime where the chemical potential is negative and |p| 3> |Ap^(p)|, 
the BCS-Leggett theory gives the lower bound of 6v~^ as 


{PpSVc) ^ — 


64 

5|fco 


E 


FL,(p) 


pZ 


+ 0 


l^ol 


= 0.44 + O 


l^ol 


(34) 


Strictly speaking, although the TMA result for coincides with Eq. fl3^ in the weak- 
coupling limit, it is still different from Eq. (1341) even in the strong-coupling limit. This is 
because the hnite value of the effective range (/cq = —30pp) causes an effective repulsive 
interaction between bound molecules in the latter limit, leading to the so-called quantum 


depletion 


Indeed, in the strong-coupling BEC regime, the last term in Eq. fl26|l . which 
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FIG. 7: (Color online) Superfluid phase transition temperatures and for various values 


of n in the cutoff function Fn{p). We take (ppU 


.3,, \-i 


= 0 . 


describes fluctuation corrections to the mean-field Green’s function Go(p, modifies the 
mean-held number equation as 


iVp 1 \ 

0 ~ A 2 -^ 




3^/7^ > 

i=x,y,z 


E ( 


(35) 


where aB,x = (374/15) |fco I and aB,y = ub.z = clb,x/3. The last term in Eq. fl3^ has 
the same form as the quantum depletion in a Bose superhuid with Nf/2 bosons, when we 
interpret 03,2 as an effective repulsive interaction between tightly bound pj-wave molecules. 
When we include this quantum depletion, the lower bound in Eq. (IHTll is improved as 

157 ^ V‘^rn\p\ 


(Pf^Vc) = —— + O 


(36) 


1357r ' \ l^ol y \ l^ol 

which agrees well with the TMA result (solid line in Fig. [ 6 ]) in the strong-coupling regime. 

Before ending this section, we comment on the the cutoff function Fn=^{p) we are using. 
As shown in Fig. [7l while the phase transition temperature is almost independent of n, 
j,px+ipy parameter. Generalizing Eq. (|3Tll to the case with an arbitrary n, 

one hnds that the lower bound (in the BGS-Leggett theory) explicitly depends on n as 

\[2, 


(pidv^) = 






im\ 


(37) 


Idvr 2n) \ 2nJ ' \ |/co| 

These n-dependences of and 6v~^ are because the factor in Ap^(p) = PxbxFn{p) 

enhances this superhuid order parameter in the high momentum region, so that physical 
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FIG. 8: (Color online) Calculated angle-resolved density of states (ARDOS) p{uj,9p). (al)(a2) 
= -12. (bl)(b2) {pIv^)~^ = -8. (cl)(c2) = -4. (dl)(d2) = 0. Upper 

and lower figures show the results at Ti"" and , respectively. In each hgure, we offset the 

results by 0.2. 

quantities in the pa;-wave superfluid phase depend on how this enhancement is suppressed 
by the cutoff function Fn{p). This implies that, in addition to the observable parameter 
set fco), one need one more experimental information about the high momentum 

regime of a real p-wave interaction, in order to unambiguously predict the phase boundary 
between the p^-wave and -|-ip^-wave superfluid phases. We briefly note that Fig. [7]shows 
that our choice (n = 3) is close to the case of discrete cutoff (which corresponds to n = cxd.) 

IV. ANGLE-RESOLVED DENSITY OF STATES AND STRONG COUPLING EF¬ 
FECTS NEAR THE PHASE BOUNDARIES. 

Figure [H] shows the angle resolved density of states p{u, 9p) (ARDOS) at (upper 
figures) and (lower figures). In Fig. [Hl(al), a dip structure is shown around w = 0. 

Since the superfluid order parameter vanishes at this is just a pseudogap originating 
from p-wave pairing fluctuations. This many-body phenomenon is non-monotonic in the 
sense that, while this pseudogap is more remarkable in Fig. IH](bl), it gradually becomes 
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obscure with further increasing the interaction strength, to eventually vanish, as shown in 
Figs. IHcl) and (dl). Figures [8^al)-(cl) also show that the pseudogap structure at is 
anisotropic in momentum space, and is the most remarkable in the direction (cos 6*p = 1). 

To simply explain this anisotropic pseudo gap phenomenon, it is convenient to employ the 
static approximation for pairing fluctuations Noting that the particle-particle scattering 
matrix F“+(q = 0, = 0) in the Px-waye channel diverges at j^, we may approximate 

the (1,1) component of the TMA self-energy in Eq. ffT4)) in the normal state near to 

O _ 

Eu{p,iUJn) 


q,tUn 


iun + L 


(38) 


where 


Apk(p) = 


-V r- 

f3 ^ ^ 


[q,lUn 


q,lVn 


VlKiP) = K.^PlF'n{P) 


(39) 


is the so-called pseudogap parameter [40|,l42|]. In obtaining Eq. fl38l) . we have only retained 
effects of the strongest px-wave pairing fluctuations near and have ignored fluctuation 
contributions from the Py-wave and p^-wave Cooper channels. Substituting Eq. fl38l) into 
the (1,1) component of the TMA Green’s function in Eq. (IT^ . one obtains 

1 


GniP,ia)n) = 


IbJr, 




^Pg(P) 


'^OJn "F 

iUJn T 


(40) 


+ + ^pgip) 

The first line in Eq. (HOj) means that the pseudogap parameter Apg(p) works as a coupling 
between the particle branch u = ^p and the hole branch u = —^p. On the viewpoint of this 
particle-hole coupling, the pseudogap may be interpreted as a result of the level repulsion 


between the particle and hole branches around a; = 0 33|, l4l|, |45|. 


The last expression in Eq. (HD]) is just the same form as the diagonal component of 
the Green’s function in the ordinary mean-held BCS theory. This coincidence immediately 
gives the BCS-type single-particle excitation spectra + l^pg(p)P) having the 

excitation gap. 


AE{ep) = 


2|fepg cos 6p I ^2 m/i — |cos 

2\p\ 


(/i > m|6pg cos6'pp), 
(/i < m|fep cos^pp). 


(41) 
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where the cutoff function Fn=z{p) has been approximated to unity. (Note that pc 3> p-p-] 


This gap is actually a pseudogap, when one correctly indue 
Cooper pairs, which is ignored in the static approximation 


es a hnite lifetime of preformed 


42|. 


Equation fITT]) indicates that, as expected, the anisotropic pseudogap phenomenon shown 
in FigJH](al)-(cl) originates from the anisotropic Px-’waye pairing fluctuations, described by 
the pseudogap parameter Apg(p) oc pi- When Ux = Uy = U^, fluctuations in all the three 
Pi-wave Cooper channels (z = x, y, z) are equally enhanced near the superfluid instability, 
so that the pseudogap parameter in Eq. (I3^ is replaced by the isotropic one. 


= 


- V 

^ / V a;,a 




q,^Vn 


P F„(.P) 


(42) 


where we have used the symmetry property, T^ j: = T^ + = The resulting pseudogap 


is isotropic in momentum space 


y,y 


Z^Z ' 


33|. 


Equation flTTl) also shows that the (pseudo)gap size AE{9p) becomes isotropic in the 
strong coupling regime where the Fermi chemical potential is negative 20|, l2l| . This is be¬ 
cause most Fermi atoms form tightly bound molecules in the strong-coupling regime, so that 
the threshold energy of single-particle excitations is simply dominated by the dissociation of 
these molecules with the binding energy Ebind — 2|/z|, as in the strong-coupling BEC regime 
of the s-wave case 50-^. 

The pseudogap parameter Apg(p) in Eq. (I39|) also explains the non-monotonic behavior 


of the pseudogap structure in terms of the interaction strength shown in Figs. [8](al)-(dl) 


33l |. Since pairing fluctuations are stronger for a stronger pairing interaction, the factor 


6pg appearing in Eq. fl39l) also becomes larger, which enhances the pseudogap parameter 
'^pg(p)- At the same, since strong pairing fluctuations are known to decrease the Fermi 
chemical potential p [20, l2l|] as shown in Fig. [3l the effective Fermi momentum defined 
by pF = becomes small. This decreases the pseudogap parameter at the effective 

2mp. As a result, while the pseudogap hrst 
becomes remarkable with increasing the interaction strength in the weak-coupling region 
because of the enhanced pairing fluctuations, it gradually shrinks when the decrease of the 
Fermi chemical potential dominantly contributes to Apg(p). In the case of Fig. [Hl(dl), one 
has p{T^"') ~ 0. Thus, the low momentum region |p| ~ 0 dominantly contributes to the 


Fermi momentum because Apg(pF) ~ 


Pi, 


density of states around a; = 0, leading to the vanishing pseudogap in this hgure 


M|. 


We briefly note that the vanishing pseudogap in the intermediate coupling regime 
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0) is quite different from the s-wave case, where the pseudogap monotoni- 
cally develops, as one passes through the BCS-BEC crossover region. This is simply because 
the contact-type s-wave pairing interaction is independent on the momentum p, so that the 
factor px is absent in the s-wave pseudogap parameter. 

At since the p^-wave superfluid order parameter Ap^{p) = bxPxFn=^{p) oc p^ 

is already present, ARDOS p{u, cos 9p = 1,0.5) in the low energy region is dominated by 
the superfluid energy gap, as shown in Figs. [8])a2)-(c2). On the other hand, such a gap 
structure is not shown in p{oj, cos9p = 0) in the weak-coupling regime (Fig. [8](a2)), because 
of the vanishing p^-wave superfluid order parameter there. However, ARDOS in the nodal 
direction (cos 6*p = 0) gradually exhibits a dip structure around a; = 0, with increasing the 
pairing interaction. (See Figs. [Hl(b2) and (c2).) Since the p^ + ipy-wave superfluid order 
parameter still vanishes at ^ this is a pseudogap induced by fluctuations in the Py- 

and p^-wave Cooper channels. Indeed, when we apply the static approximation to the region 
near pP^'^^Py^ the (1,1) component of the single-particle Green’s function with p = {0,py,pz) 
is reduced to Eq. fjTnil where the pseudogap parameter is replaced by 


^lg{^,Py,Pz) = 


i=y,z 


/S 




fqr, iiy„ 


q.ll'n 


plFlip) = X] ypg^PiFLzip)- 


(43) 


i=y,z 


Although Eq. fjTHll is similar to Eq. ([39]), the former involves effects of px-wave superfluid 
order parameter Ap^{p). Since gapless Fermi excitations only remain along the line node of 
the Px-wave superfluid order parameter, pairing fluctuations described by 
in Eq. fH3|) are weaker than pairing fluctuations described by &pg(T = Tp^) in Eq. (I39|) . 
This explains why the pseudogap appearing in ARDOS p{u,cos9p = 0) at (Figs. 

IH](b2) and (c2)) is less remarkable, compared to the dip structure in p{u;, cos 9p = 1) at Tp^ 
(Figs. IH](al)-(cl)). 

The reason for the vanishing superfluid gap and pseudogap gap in Fig. IH](d2) is the same 
as that in the case of Fig. [8](dl). That is, at this interaction strength, the chemical potential 
is very small (/i(T(f“’’''*^^) = O.OSep), so that the low-momentum region (|p| ~ 0) dominantly 
contributes to ARDOS p{oj, 9p) around oj = 0. Thus, the Px-wave superfluid order parameter 
Ap^ip) oc Pxi as well as effects of the pseudogap parameter Apg(0,Py,p^) cc. p^ + pi in Eq. 
fH3l) . do not almost affect ARDOS around a; = 0 in this case. 

Figures [9](al) and (a2) show that the pseudogap in p{oj, cos 9p = 1) in the normal state 
continuously changes to the superfluid gap, as one passes through the superfluid instability 
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FIG. 9: (Color online) Calculated angle-resolved density of states (ARDOS) p{uj, 9p) at various 
temperatures, when ((ppUa:)”^, (pp<5n)“^) = (—8,0.3). (al)(a2) cosOp = 1. (bl)(b2) cos0p = 0.5. 
(cl)(c2) cos Op = 0. The upper and lower figures show the results in the normal state, and 
in the Px-wave superfluid state, respectively. Evaluating the pseudogap temperature T*{0p) as 
the temperature below which a dip structure appears in Op) around a; = 0, one obtains 
T*{cos0p = 1) = 0.13Tp, T*{cos0p = 0.5) = O.llTp, and T*{cos0p = 0) = O.lTp. In each figure, 
we offset the results by 0.3. 


at The same phenomenon is also shown when cos 6p = 0.5, as shown in Figs. [9](bl) 
and (b2). On the other hand, since the px-wave superfluid order parameter vanishes when 
cos 9p = 0, Figs, in] (cl) and (c2) show how the pseudogap in the nodal direction continues 
developing in the px-^ave superfluid phase, to be the most remarkable at . 

When we introduce the characteristic temperature T*{6p) as the temperature below which 
a dip structure appears in p{oj, cos 9p), we obtain Fig. [TOl Because Ap^(p) = 0 in the 
nodal direction (cos^p = 0), we may regard T* {cos 9p = 0) as the pseudogap temperature 


33l. l4ll. l45l | in this momentum direction, below which strong pairing fluctuations induce a 


pseudogap in ARDOS p{u, cos 9p = 0). In this case, one may call the region surrounded by 
T* {cos 9p = 0) and the pseudogap regime. 

In the case of cos 9p ^ 0, T*{9p) also has the meaning of the pseudogap temperature, when 
T*{9p) > T'P^. On the other hand, T*{cos9p = 1, 0.5) is lower than around (plua,)”^ = 0 
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FIG. 10: (Color online) Characteristic temperature T*{9p) below which a dip structure appears 
in ARDOS p{io,9p). We take {ppSv)~^ = 0.3. The dashed-dotted line shows the temperature at 
which the Fermi chemical potential p vanishes. The chemical potential p is negative in the right 
side of this line, so that this strong-coupling regime may be regarded as a gas of two-body bound 
molecules with the binding energy Sbind ~ 2|/r| 33|, rather than a gas of Fermi atoms. 


in FigJTOl which means that the superfluid gap does not appear in ARDOS when T*{6p) < 
T < As mentioned previously, since \p\ <C in this intermediate coupling regime, 

single-particle excitations around p = 0 dominantly contribute to ARDOS around u = 
0. Because of this, a small superfluid excitation gap by a small p^-wave superfluid order 
parameter, Ap^(p) ~ bx^/2mp ~ 0, around p = 0 is easily smeared out by strong pairing 
fluctuations existing in this regime even below Since this strong-coupling effect is 

gradually suppressed below ARDOS starts to exhibit a superfluid gap structure below 
T*{6p) (See Fig. [TTf a).). to approach the BCS-type superfluid density of states shown in 
Fig. fTlT b). Thus, T*{6p) in this regime may be regarded as the characteristic temperature, 
below which the p^-wave superfluid order overwhelms pairing fluctuations. 

In Fi g. [Ihl we also plot the temperature at which the Fermi chemical potential p changes 
its sign As mentioned previously, in the strong-coupling regime where p < 0, since 
the system is dominated by tightly bound molecules, the p-wave character of Cooper pairs 
is less important. In the normal state near this fact gives the isotropic pseudogap 
size AE{6p) in Eq. fjTTl) . In the px-wave superfluid phase, the Bogoliubov single particle 
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FIG. 11: (Color online) (a) Angle-resolved density of states p{u,cos9p = 1) in the px-wave super¬ 
fluid phase. We take ~ (0,0.3). Each line is offset by 0.03. (b) Mean-field 

result at T = , which is obtained by ignoring the self-energy correction S(p,zu;„) in Eq. 

(fT^ in calculating ARDOS. 

excitation spectrum in the strong coupling regime, 

Ep = \J {£p + \x\y + \^pApW ) ( 44 ) 

also has the isotropic energy gap 2|/x|, reflecting that the threshold energy of Fermi ex¬ 
citations is simply dominated by the binding energy (i^bind ~ 2|/r|) of a two-body bound 
molecule. Thus, the p-wave anisotropy is not important in the right side of this line in Fig. 
m as far as we consider low-energy single-particle excitations. 

V. SUMMARY 

To summarize, we have discussed strong-coupling properties of a one-component super¬ 
fluid Fermi gas with a uniaxially anisotropic p-wave pairing interaction (Ux > Uy = Uz)- 
Including p-wave pairing fluctuations within a T-matrix approximation, we determined the 
two superfluid phase transition temperatures which gives the phase boundary between 
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the normal state and the Px-wave superfluid state, and (^< T|'^), which gives the 

phase boundary between the Pa,-wave and + ipy-wave superfluid states. 

We examined single-particle excitations near as well as near ^ In the normal 

state near we showed that strong pairing fluctuations in the Pa;-wave Cooper channel 
induce an anisotropic pseudogap phenomenon where a pseudogap structure in the angle- 
resolved density of states (ARDOS) is the most remarkable in the px direction. We also 
showed that this pseudogap continuously changes to the p^-wave superfluid gap below . 
On the other hand, the pseudogap was found to continue developing below in the nodal 
direction (_L px) of the Px-wave superfluid order parameter, to be the most remarkable at 
rjnVx+tpy ^ Since pairing fluctuations are simply suppressed in an isotropic s-wave superfluid 
state, this phenomenon is characteristic of a p-wave Fermi superfluid with a nodal superfluid 
order parameter and with plural superfluid phases. To characterize the anisotropic pseudo¬ 
gap phenomenon, we determined the characteristic temperature T*{6p), below which a dip 
structure appears in ARDOS. 

In this paper, we have considered the normal state, as well as the Px-wave superfluid phase. 
To obtain the complete understanding of a p-wave superfluid Fermi gas with a uniaxially 
anisotropic p-wave interaction, we need to also examine the px + ipy-wave superfluid phase 
below In addition, for simplicity, we employed the BCS Hamiltonian with a p-wave 

pairing interaction, which implicitly assumes a broad Feshbach resonance. In this regard, all 
current experiments are using a narrow p-wave Feshbach resonance, so that it is an important 
problem to clarify how the resonance width affects strong-coupling properties of a p-wave 
superfluid Fermi gas. Since a p-wave superfluid state is known to be sensitive to spatial 
inhomogeneity, inclusion of a realistic harmonic trap also remains as our future problem. 

The realization of a p-wave superfluid Fermi gas is an exciting challenge, in order to 
qualitatively go beyond the current stage of cold Fermi gas physics that the s-wave Fermi 
superfluid has only been realized. Since the anisotropic pairing is a crucial key in a p-wave 
Fermi superfluid, our results would contribute to understanding how this character affects 
many-body properties of an ultracold Fermi gas, especially in the superfluid phase. 
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